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Abstract: The enhancement of pions and kaons observed at small transverse momenta 
in ultra-relativistic heavy-ion collisions may at least partly reflect critical scattering expected 
1. . to occur in the neighborhood of a second order phase transition. Kinetic equations in the 
relaxation time approximation are proposed for the time evolution of the quark distribution 
function into that of the pions. Relaxation times for thermalization and hadronization 
processes are functions of momenta and approach zero in the limit p — > 0, a consequence of 



X 



criticality at the phase transition. Data can be reproduced for suitably chosen parameters. 
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Introduction 

The low-pT enhancement observed in ultra-relativistic heavy-ion collisions is one of the 
unexpected results in this domain of physics (c.f. the review PJ): At small values of the 



transverse mass, mx = Jm 2 + pj,, the observed distributions of pions and kaons, dN/dmT, 
considerably exceed the thermal distribution exp(— mr/T) @. Various explanations have 
been proposed: e.g. existence of the transverse flow [^ ffl, creation of small plasma droplets 
0, decays of resonances [||, formation of the pion system out of chemical equilibrium |7|, ||, 
or the medium modification of the pion dispersion relation ||. In this paper we propose yet 
another explanation, namely that the low-pr enhancement reflects critical scattering. 

Critical scattering or critical opalescence is a rather general phenomenon observed in 
the neighborhood of a second order phase transition ||L0 |: The cross section d 3 a/dc[ for the 



scattering of light, X-rays or neutrons on a medium is given by the expression 

d 3 a 



dq 



oc 



f d 3 r /d 3 rV q(r - r ' ) (n(r)n(r')) T , (1) 



where q is the momentum transfer and (n(r )n(r' ))? is the density-density correlation func- 
tion for the particles in the medium at the temperature T. In the neighborhood of a phase 
transition a long-range ordering develops and the cross section for small values of q has the 
form 

— - (2) 

dq K P(T) + q 2 ' l ' 

where 

k\T) oc | 1 - L r (3) 

J-c 

is the inverse correlation length with 7 being the critical exponent. 

This kind of experiment, where an external test particle is scattered on the system in its 
critical state, cannot be performed for a quark-gluon plasma which only lives for 10~ 23 s. 
However, also particles which are created inside the plasma should experience critical scat- 
tering and their observed momentum distribution may bear signal of this phenomenon. In 
the present paper we discuss one possible realization of such a scenario. 



In explicit calculations [11, 12 1, using the Nambu - Jona-Lasinio (NJL) model for a 
description of the quark-meson plasma (no gluons), it has been shown that singularities 
really occur in the cross sections at small center-of-mass energies y/s of colliding particles, 



when the phase transition is approached. In particular, the integrated elastic quark-antiquark 
cross-section a qq ^ qq (s, T) diverges like s _1 for T —>■ T c . A singularity is also found for the 
hadronization cross-section qq — > im. In both cases the singularity arises because the quark 
condensate (qq)r, which is the order parameter of the chiral phase transition, goes to zero 
at the critical temperature. 

May the observed low-pT enhancement originate from the low y/s singularity in the cross 
sections? We think this is so and demonstrate it using a toy model. Kinetic equations are 
proposed for the time evolution of a spatially homogeneous quark plasma (described by the 
single particle quark distribution function f q (p,t), p = |p|) into a pion gas (with a corre- 
sponding distribution function f^ip, t)). Criticality of the system enters via the singularities 
in the thermalization and hadronization cross sections. 

In order to keep the model as transparent as possible we work in the relaxation time 
approximation, thus reducing the kinetic equations to a coupled system of first order integro- 
differential equations for f q (p,t) and f-nip.t). Contrary to the usual approaches, the relax- 
ation times r(p, t) depend on momenta via the energy dependence of the cross sections. The 
singularities in the cross sections due to criticality of the medium translate themselves into 
singularities of the inverse relaxation times l/r(p,t) in the momentum p. 

Kinetic description of the hadronization of the quark plasma 

The above ideas are translated into the following set of kinetic equations 



dfq(p,t) fq(p,t) - fq,th(p,t) _ fq(j>,t) fir(p,t) 

dt T th (p,t) r had (p,t) r dec (p,t) 



dfnfat) fn(p,t)~ fq(p,t) _ f n (p,t) _ fn(p,t) ,_ 

dt r th (p,t) T had {p,t) r dec (p,t) T em (p) ' 

dfem(P,t) fn(P,t) 



dt r em (p) 



(6) 



Eq. (U) determines the time evolution of the quark distribution function f q (p,t) (it de- 
scribes both quarks and antiquarks, i.e. f g (p,t) = f qU arks(p,t) + fanUquarks(p,t))- The first 
term on the r.h.s. of Eq. (|j) is the collision term written in the relaxation time approxima- 
tion; it is responsible for the thermalization of quarks since the distribution function f q (p,t) 



is always attracted to the thermal one f q ,th(p,t). The second (third) term describes the loss 
(gain) of quarks due to the hadronization (deconfmement) process where we limit ourselves 
to the reaction qq — > tttt (qq <— Tin). Of the possible reaction channels qq — > nn (n > 2), the 
case n = 2 may be the dominant one for phase space reasons at moderate energies in the qq 
system, and we assume that this is so. 

In Eq. (J5|), for the time evolution of the distribution functions of pions, the first three 
terms on the r.h.s. describe: the thermalization of pions, appearance of pions due to 
hadronization of quarks, and the two-pion reaction into quark-antiquark pairs, respectively. 
The last (fourth) term accounts for the emission of pions from the plasma, which in a realistic 
system proceeds via the surface, while we simulate it via a homogeneously distributed sink. 
The emission of preequilibrium pions is crucial in order to observe deviations from thermal 
equilibrium. 

Eqs. (fD - (H) can be solved for any initial conditions (i.e. assuming some particular form 
of the distribution functions at the initial time t — 0) provided the thermal distribution 
functions f q ,th(p,t) and are known at all times. Since the first terms on the r.h.s. of Eq. 
@ and (|D represent the collision terms (in the relaxation time approximation) they must 
obey the symmetries leading to particle and energy conservation. This gives the following 
constraints for f Qt th(p,t) and at each time t 



d 3 P fi(p,t) ~ fi,th(Pit) 



(2tt)3 Tth (p,t) 



d3 p nrz 2 fi(p>t) - fi,th(p,t) 



(7) 



p 1 + m. 



(2vr)3V- . TthM 



(here i = q or 7r). Eqs. (|7|) and (|8|) determine the temperature Tj(t) and the chemical 
potential /!«(£) appearing in the thermal distributions. For simplicity we assume the low- 
density high-temperature form of these functions, namely a Boltzmann distribution 



fi,th(p,t) =5 , i exp 



yV + ml - IM(t) 
Ti(t) 



(9) 



where gi are the degeneracy factors: g q = 24 (quarks and antiquarks having two different 
spin projections, 3 colors and 2 flavors) and g n = 3 (three different values of the isospin). 



Relaxation times and cross sections 

The analysis of the exact collision term in the Boltzmann kinetic equation (analogous to 
that in |13| ) leads to an expression for the average time for hadronization of two quarks into 
two pions 

1 1 f d 3 Pl 



/d p-\ 
— — nr— 2 /,(Pi, t)F„(8)**^{s), (io) 



where a q q_> nn (s) denotes the total cross section for this process. The relativistic flux factor 



of incoming quarks is F qq (s) = |v/s(s — 4m 2 ,), with y/s being the center-of-mass energy of 
the quarks with momenta p and pi, respectively. The expression giving the deconfinement 
relaxation time, Td ec {p,t), has the form analogous to Eq. (|10|). 

The relaxation time for the thermalization of quarks can be written in the form 



--f q (pi,t)F qq (s)a qq ^ qq (s) 



Tth(p,t) 2yjp> + m 2 q J (27i) 3 ^pj + m 2 q 

I r d V\ 

+ o / 2 , 2 / /n Vt / 2 =A(Pl»0^«pr(a)^-^(g)» (11) 

2Wj9 2 + m^ J (27r) 3 i/p| + m^ 



where the two terms on the r.h.s. of Eq. ( |TT1) appear due to the quark-quark and quark- 
pion elastic collisions respectively. The formula for the thermalization time of pions can be 
obtained formally from Eq. ( pTTj) by the exchange of the indices q and it. 

For the sake of simplicity, in what follows we shall treat the quarks and pions as massless 
particles. We shall also assume that the cross sections: cr qq ^ n7T (s),a qq ^ qq (s) : a nn ^ 7T7r (s) and 
0"<j7r^<j7r(<s) are given by the generic expression 



a(s) = Co 



1 + ^ 
s 



(12) 



where <To represents a constant contribution to the cross sections, and the appearance of the 
singular term sq/s accounts for the critical phenomena. The hadronization and deconfine- 
ment cross sections, a qq ^ 7T7T (s) and cr n7T ^ qq (s) , are related to each other by the principle of 
detailed balance 



g q a m ^ im {.s)=gi t <7i m -*qq(s). (13) 

This relation guarantees that in the absence of emission (r em — > oo), the kinetic equations 
lead to the chemical equilibrium riq/n^ = g q j ' g v . 

Since the relaxation time approximation requires that the system is always close to ther- 
mal equilibrium, the distribution functions appearing in Eqs. (JT^ - [TTj) can be replaced by 
the thermal ones and, in this way, we find the approximate form for the hadronization time 



a n q (t) 



T ha d(p,t) 
where 



1+ S ° 



ApT q (t) 



(14) 



d 3 p 

(2^ 



«*(*) = ItA-Mp^)- (15) 



Eq. ( pf ) can be used to find the characteristic momentum scale, Sp(t) = s /AT q (t), below 
which the hadronization time significantly deviates from a constant. Due to the faster 
hadronization of low-energy quarks their temperature T q (t) increases and, correspondingly, 
5p(t) decreases in time. Therefore, in the following we shall use the time average of Sp(t) , i. e. 
the quantity 5p = (5p(t)), in order to make the estimate of the momentum range, < p < Sp, 
where the enhancement in pion production is expected. 

In our model the pions are emitted from the whole volume of the interacting system. 
The most natural physical assumption for T em (p) would be, that it is inversely proportional 
to the velocity of the pion, which for massless particles equals c for all momenta. For this 
reason we choose r em to be a constant. In order to study the sensitivity of our results on the 
emission rate we shall introduce the ratio TZ defined as 

1 1 

— = n — ( r~m- ( 16 ) 

Tern T t h{P ~> OO, t = 0) 

The magnitude of TZ indicates how much faster the evaporation is, in comparison to the rate 
of the thermalization processes inside the system. 



Results 

In order to reduce the number of free parameters in our calculation we have set all the 
cross sections: <J qq ^Tnr(s) , cr qq ^ qq (s) , cr qn ^ qn (s) and ov^^^s) to be equal (within a factor of 
3 such a result is supported by calculations within the NJL model). Then, we are left with 
three parameters: ctq, Sq and 1Z. One can notice that <To sets the overall time-scale and is 
unimportant if we are interested in the final (t — > oo) pion distributions. 

We start solving our kinetic equations by assuming that the initial quark distribution 
is a thermal one and that there are no pions in the system. After integrating Eqs. M) - 
(|6|) till the time when all quarks are hadronized and all pions are emitted, one obtains the 
distribution function f em (p, t — > oo) of the observed pions for a set of parameters s and 1Z. 

In Fig. 1 we show our results for the case when the initial temperature of the quarks 
T q (t = 0) = 140 MeV. The dashed lines represent the initial quark distribution functions, 
whereas the solid ones represent the distribution functions of the emitted (observed) pions for 
four different choices of the parameters: (a) 1Z = 33, s = 0.48 GeV 2 ; (b) 1Z = 11, s = 0.48 
GeV 2 ; (c) 1Z = 33, s = 0.16 GeV 2 ; and (d) 1Z = 11, s = 0.16 GeV 2 . The corresponding 
values of the momentum scale Sp are: 210, 260, 90 and 95 MeV, respectively. 

In all the considered cases one can clearly see the enhancement in the pion production 
for p < 5p. During the evolution of our system the low energy quarks hadronize more 
effectively because of the p-dependence of Thud and, correspondingly, the low energy pions 
are preferably produced. Our results also show that the effect depends strongly on the pion 
emission rate: for larger values of TZ the effect is more significant. Such a dependence can be 
easily understood since the thermalization leads always to the distribution functions of the 
form (|^), and in the limiting case 1Z — > no effect could be seen. The last fact indicates also 
that our mechanism for obtaining the low-pr enhancement is basically a non-equilibrium one: 
The excess can be observed only if the pions are emitted from a non-equilibrium system. 

In Fig. 2 we plot the ratio r(p) obtained by normalization of the distribution function 
of the observed pions to the Boltzmann distribution. In this case the latter is obtained by 
fitting the exponent function to f em (p,t — > oo) in the region 0.5 GeV < p < 1 GeV. The 
values of the parameters are the same as those assumed in Fig. la. Our result is shown 
together with the data 0] (14.6 A GeV Si + Pb -» vr", rapidity interval y = 3.4 - 3.6). 
Because we did the calculations in the special case m n = 0, we treat the data as if they were 
obtained also for massless pions and identify the quantity itlt — m n (appearing in [14]]) with 



the transverse momentum px- Moreover, using the relation px = \/2p/3 (valid for isotropic 



systems) we find m T — m^ = J2p/3. 

Fig. 2 shows that our model is able to explain a characteristic shape of the observed 
enhancement, if the two parameters (sq and TV) are adjusted. Nevertheless, we are of the 
opinion that the low-pr enhancement in pion production is caused by the superposition 
of various mechanisms. In Ref. fLjJ it is argued that the decay of A resonances can be 



responsible for the excess of low-p^ pions, although in different rapidity intervals different 
ratios of pions from A decay to direct pions must be assumed to fit well the data. This 
uncertainity leaves room for other effects. Also the origin of the enhancement in kaon 
production, starting at much smaller values of px, is still not clear. Therefore, we think 
that our mechanism for obtaining the low-pT enhancement is probably only one of a few 
contributions leading to this effect. 
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Figure Caption 

Fig. 1 The initial quark distribution functions (dashed lines) and the distribution func- 
tions of the emitted pions (solid lines) for four different choices of the parameters: (a) 
K = 33, s = 0.48 GeV 2 ; (b) K = ll,s = 0.48 GeV 2 ; (c) K = 33, s = 0.16 GeV 2 ; and (d) 
K= 11, s = 0.16 GeV 2 . 

Fig. 2 The distribution function of the observed pions normalized to the Boltzmann 



distribution. The parameters as in Fig. la. The diamonds show the data |T4| 
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Fig. 2 
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